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FORWARD 


In  energy  systems,  the  problems  associated  with  high  tempera¬ 
ture  corrosion  especially  in  mixed  oxidants,  e.g.  oxygen  containing 
sulfur  dioxide  and  sulfur  trioxide,  have  in  many  instances  limited 
the  use  of  certain  materials.  As  an  example  of  the  drastic  changes 
in  reaction  rate,  the  oxidation  rate  of  nickel  in  oxygen  plus  some 
sulfur  containing  gas  may  exceed  the  rate  in  pure  oxygen  by  as  much 
as  10®!  Accordingly  research  to  determine  diffusivity  of  sulfur 
through  oxides  was  undertaken  along  with  a  parallel  study  on  the 
kinetics  of  corrosion  of  nickel  in  S02“02  gas  mixtures. 


---  ■  daaaiiflrtMiiiEiflttiMiM 
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STATEMENT  OF  THE  PROBLEM 


Many  high  temperature  corrosion  processes  are  determined  by 
the  rate  at  which  an  impurity  anion  (not  oxygen)  is  transported 
through  a  corrosion  product  layer.  An  example  is  the  high  temp¬ 
erature  corrosion  of  metals  and  alloys  in  mixed  oxidants  such 
as  S02“02  §as  mixtures.  Accordingly,  research  was  initiated  in 
(1)  Corrosion  of  nickel  in  S02"02  gas  mixtures  at  800  .  (2)  Defect 
structure  of  calcia  and  of  yittria  stabilized  zirconia;  and  dif¬ 
fusion  of  radiotracer  sulfur-35  and  tellurium  127m  in  calcia 
stabilized  zirconia.  The  tracer  studies  bear  directly  on  the 
problem  of  sulfur  penetration  into  chemical  sensor  materials  and 
also  the  penetration  of  sulfur  through  thermal  barrier  layers  of 
zirconia  on  turbine  blade  alloys. 
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ABSTRACT 


This  study  was  divided  into  three  interrelated  sections. 

These  were  (1)  corrosion  of  nickel,  in  SO2“02“Ar  gas  mixtures 
at  800°c ;  (2)  defect  structure  of  calcia  and  of  yttria  stabilized 
zirconia;  and  (3)  diffusion  of  radio-tracer  sulfur-35  and  of 
tellurium  127m  in  calcia  stabilized  zirconia.  The  corrosion 
study  resulted  in  a  confirmation  of  linear  kinetics  for  cor¬ 
rosion  of  nickel  in  S02"Argon  gas  mixtures  at  800°C  and  an  agree¬ 
ment  of  a  self-consistent  mechanism  involving  the  migration  of 
SC>2  as  a  gas  molecule  through  the  outer  NiO  layer  and  dissociation 
at  the  inner  sulfide  layer  as  a  rate  determining  step.  The  cor¬ 
rosion  kinetics  in  20%02~variable  S02"balance  argon  and  in  30% 

SO2  variable  02~balance  argon  were  extremely  rapid.  These 

data  did  not  show  that  interconnected  stringers  of  sulfide  through 

the  outer  NiO  layer  was  responsible  for  the  observed  kinetics. 

The  ionic  conductivity  and  diffusion  in  calcia  stabilized 
zirconia,  in  yttria  stabilized  zirconia  and  related  structures, 
were  developed  using  a  theoretical  model  which  when  fitted  to  the 
existing  expermental  data  yielded  good  agreement.  Briefly,  the 
increase  in  conductivity  with  dopant  concentration  followed  by  a 
maximum  and  a  subsequent  decrease  in  conductivity  wa s  explained  by 
considering  jumps  of  oxygen  via  oxygen  vacancies  one-fold  and  two¬ 
fold  coordinated  calcium  (or  other  aliovalent  dopants)  in  the 
zirconia. 

Radiotracer  sulfur-35  or  tellurium-127m  was  diffused  into 
15  m/o  Ca0-Zr02  (CSZ).  A  multimode  diffusion  process  was  observed. 
The  two  bulk  diffusion  processes  were  related  to  jumps  via  anion 
vacancies  in  the  CSZ.  Surprisingly,  there  was  little  or  no  cor¬ 
relation  between  size  or  mass  in  these  studies,  i.e.  sulfur  and 
tellurium  both  diffused  at  approximately  the  same  rate.  The  dif- 
fusivities  were  about  four  orders  of  magnitude  greater  than  sul¬ 
fur  in  NiO  in  agreement  with  the  large  vacancy  concentrations  in 
CSZ.  Furthermore,  both  sulfur  and  tellurium  diffused  more  slowly 
than  oxygen  but  faster  than  either  Ca  or  Zr  in  CSZ. 
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L.  Corrosion  of  nickel  in  S02"02-Ar  Gas  mixtures  at  800°C. 

Polycrystalline  nickel  samples  (99.6  +  w/o)  approximately  25 
by  17  and  0.5mm  were  polished  mechanically  on  SiC  paper  from  240 
grit  to  600  grit.  Some  samples  were  pre-annealed  in  one  atmosphere 
of  hydrogen  at  about  800°C  prior  to  the  experiment.  In  addition, 
a  few  experiments  were  carried  out  on  single  crystals  of  nickel 
(99.99  +  w/o).  Results  were  the  same  as  for  polycrystalline 
samples.  Mixtures  of  S02“02~Ar  were  obtained  using  metered  cap¬ 
illary  flowmeters.  The  gas  mixtures  are  indicated  on  the  stability 
diagram  in  Fig.  1.  Kinetics  were  determined  using  an  Ainsworth 
automatic  recording  balance.  Preliminary  experiments  were  per¬ 
formed  on  nickel  in  pure  oxygen  to  test  for  self  consistent  data. 

See  Fig.  2. 

A.  Corrosion  in  S02“Argon  gas  mixtures. 

At  800°C,  the  kinetics  in  S02"Ar  gas  mixtures  were  linear. 

See  Fig.  3.  These  data  are  in  good  agreement  with  the 
£ata  of  Nakai  and  Takazawa  (1)  and  those  of  Kofstad  and 
Akesson  (2).  See  Fig.  4.  According  to  (2),  the  rate 
determining  step  is 

(1)  S02  (g)  2  SO2  (ads) 

(2)  SO^  (ads)  Reaction  products. 

Our  data  are  in  accord  with  their  suggestion.  A  suggested 
mechanism  along  with  calculated  activities  and  partial 
pressures  are  shown  in  Fig.  5.  Markers  were  always  found 
at  the  sulfide-oxide  interface  corrosion. 

B.  Corrosion  in  207.  02  “  variable  S02  -balance  Ar.  at  800°C. 

These  kinetics  exhibited  protective  behavior  at  ^20%  S02 
and  were  exceedingly  rapid  at  s  307.  S02-  See  Fig.  6. 

The  drastic  change  in  kinetics  suggest  a  change  in  mechan¬ 
ism.  We  tested  the  resistance  of  the  scale  of  nickel 
sulfide  containing  some  sulfide.  In  no  case  did  we  find 
evidence  by  electrical  measurements  nor  by  optical  and  SEM 
metallographic  examination  of  stringers  of  sulfide  forming 
a  continuous  path  through  the  NiO  scale  as  was  reported  by 
Luthra  and  Worrell  (3,4)  for  similar  corrosion  at  603°C. 

In  both  the  cases  of  S02“Ar  and  20"  02~variable  S02”balance  Ar, 
both  the  inner  layer  of  sulfide  and  the  outer  layer  of  NiO 
(containing  sulfur  but  no  continuous  stringers)  increase 
in  thickness  with  time.  See  Fig.  7.  We  are  currently 
studying  the  conductivity  of  Ni3S2  “  NiO  two  phase  mixtures 
contribute  to  an  understanding  of  these  phenomena.  (5). 


C.  Corrosion  in  307,  SC>2“variabLe  02"balance  Ar. 

Initially  as  the  oxygen  content  is  increased  (27,,  57,,  10%) 
the  kinetics  decrease.  However,  for  s  20%,  the  kinetics 
increase  drastically.  See  Fig.  8.  At  the  higher  values 
of  P02>NiSO4  formed  on  the  outer  surface  of  the  NiO  layer 

in  accord  with  the  stability  diagram.  See  Fig.  1.  Moreover, 
the  markers  were  found  at  the  inner  sulfide  layer-outer 
nickel  oxide  layer  interface  as  in  in  the  study  mentioned 
above.  However,  at  lower  oxygen  pressures,  the  markers 
were  found  in  the  middle  of  the  outer  oxide  layer. 

D.  Effect  of  pre-oxidation  on  corrosion  mechanism. 

Clearly  the  transport  of  SO2  or  sulfur  in  some  form  through 
the  corrosion  product  affects  the  kinetics  dramatically. 
Accordingly,  nickel  coupons  were  oxidized  in  pure  oxygen 
and  then  exposed  to  sulfur  containing  gases.  See  Fig.  9. 
There  is  a  relation  between  the  thickness  of  pre-formed 
NiO  and  the  incubation  period  or  onset  of  rapid  kinetics 
after  admitting  SO2  gases.  However,  the  relation  is  not 
incubation  period  proportional  to  the  square  of  the  pre¬ 
formed  oxide  thickness  as  one  would  expect. 


2.  Defect  Structure,  Ionic  Conductivity  and  Diffusion  in  Calcia-, 

in  Yttria-stabilized  Zirconia  and  in  related  Fluorite  structures. 


Prior  to  studying  the  diffusion  of  radio  tracer  S-35  and  Te-127m 
in  Ca0-Zr02,  a  detailed  theoretical  analysis  was  performed  on  the 
existing  data  for  ionic  conductivity  of  calcia-,  or  yttria- 
stabilized  zirconia  and  other  related  fluorites.  In  the  case 
of  the  stabilized  zirconias ,  the  isothermal  ionic  conductivity 
first  increases  with  added  aliovalent  dopant,  attains  a  maximum 
and  then  decreases.  See  Appendix,  especially  Fig.  1.  While  it  is 
well-known  that  oxygen  migrates  via  vacancies  in  these  materials, 
the  occurrence  of  a  maximum  and  subsequent  decrease  was  not  well 
explained.  When  the  concentration  of  dopant  attains  12-13m/o 
(Ca0-Zr02)  and  8-9  m/o  (Y  203"Zr02 ) ,  the  anion  vacancies  are  not 
in  dilute  solution.  The  vacancies  and  dopant  atoms  form  dopant- 
defect  associates.  The  ionic  donduction  was  assumed  to  occur  via 
a  multimode  mechanism  through  and  between  Ca^' -one  fold  and  Car'  - 
two  fold  coordinated  anion  vacancy  sites.  This  model  was  succes¬ 
sfully  applied  to  calcia  stabilized  zirconia  (See  Appendix  1). 

A  similar  analysis  has  been  applied  to  yttria  stabilized  zirconia. 
Results  of  the  theory  and  experimental  data  are  shown  in  Fig.  10. 

Furthermore,  the  theoretical  model  is  also  applicable  to  various 
other  oxide  electrolyte  systems  with  the  fluorite  structure,  such 
as  the  rare  oxides  (Me203)  -  stabilized  zirconia  with  general 
formula  Me^x  Zr^^  °2-2x.  ^•Me  =  Ta»  Dy»  Sc»  Yb>  Sm»  Ho»  In 

and  also  CaO  -  and  Me20^  -  doped  CeO ^  and  ThO^  systems.  A  man¬ 
uscript  is  in  preparation  describing  these  results. 
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3.  Radiotracer  Diffusion  of  S-35  and  Te-127m  in  Calcia-stabilized 
Zirconia  and  Yttria-stabilized  Zirconia. 


Polycrystalline  samples  of  calcia  stabilized  zirconia,  (9,  13,  15 
and  18  m/o)  were  obtained  from  the  General  Electric  Company  courtesy 
D.  Chatterji  and  K.  Lay.  Polycrystalline  samples  of  CSZ  of  15  m/o 
calcia  were  obtained  from  the  Zircoa  Corporation.  Yttria-stabilized 
zirconia  polycrystalline  samples  were  obtained  from  Viking  Chemical 
courtesy  of  Dr.  M.  Kleitz  of  Grenoble  and  single  crystalline  YSZ 
(9.6  m/o)  was  obtained  from  Dr.  Joseph  Wenckus ,  President  of  Ceres 
Corporation.  In  some  experiments  the  samples  were  used  as-received 
while  in  others  the  samples  were  given  a  high  temperature  anneal 
prior  to  the  diffusion  anneal.  The  experimental  procedure  was  almost 
the  same  as  that  for  sulfur  diffusion  in  single  crystalline  NiO. 

[(Ref.  12  ).  Supported  by  ARO  under  DAAG-29-76-G-0196  and  DAH  C04- 
74-G-0130].  Cylindrical  pellets,  about  11  mm  in  diameter  by  1.7  mm 
thick  were  used  as  received  or  pre-annealed  in  1  atm  O2.  The  samples 
were  polished  on  two  flat  faces  using  SiC  paper  and  diamond  paste  of 
1  micron.  The  ionic  conductivity  of  these  samples  was  first  checked 
to  ascertain  that  the  material  behaved  in  accord  with  previous  studies. 

The  radiotracer  (S-35  in  benzene  or  Te-127m  in  4m  HCl)  was  de¬ 
posited  in  a  Pt  crucible  by  slowly  evaporating  the  solvents.  The 
Pt  crucible  and  the  sample  were  encapsulated  in  a  quartz  tube  and 
the  tube  given  diffusion  anneal  (850°C  to  1250°C  for  0.5  to  300  hrs. 
depending  on  the  temperature).  The  tubes  were  quenched,  samples 
removed  and  washed  in  acetone  or  toluene.  The  cylindrical  sides  were 
ground  away  ~  5 J  DT'  to  avoid  side  effects. 

The  tracer  penetration  profiles  were  obtained  by  counting  the 
weak  betas  and  serially  sectioning  the  samples.  The  thicknesses 
removed  were  obtained  by  weight  change.  A  high  surface  activity 
(105  _  I06  Cpm)  was  used  so  corrections  for  coincidence  loss  were 
made.  No  corrections  for  self  absorption  were  made  because  of  the 
deep  penetrations  and  weak  beta  energies. 

The  penetration  profiles  for  S-35  and  Te-127m  in  15  m/o  CaO-ZrO^ 
generally  consisted  of  four  regions:  These  were: 

I.  Near  surface  region  (ax  <5o  at  1000°C)  where  a  sharp  drop 
in  activity  is  observed. 

II.  A  first  intermediate  region  which  corresponds  to  the  diffusion 
process  with  the  slower  diffusivity  and  higher  solubility. 

III.  A  second  intermediate  range  which  corresponds  to  the  diffusion 
process  with  high  diffusivity  and  low  solubility. 

IV.  A  short-circuit  diffusion  and/or  almost  constant  activity  region 
which  is  characterized  by  a  long  "diffusion  tail"  with  much 
scatter  in  activity  for  a  given  penetration  distance. 
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Fig.  10  shows  a  typica1  penetration  profile  of  S-35  and  Te-127m  in 
15  m/o  Ca0-Zr02«  Region  I  may  be  related  to  a  surface  process.  No 
analysis  of  this  narrow  region  was  performed.  Regions  II  and  III 
appear  to  be  bulk  diffusion  processes  and  region  IV  a  short  circuit 
process  or  processes.  The  data  for  Regions  II  and  III  were  analyzed 
analogously  to  similar  data  for  NiO  (12)  and  an  Arrhenius  plot  is  shown 
in  Fig.  11. 

Because  we  used  polycrystalline  samples,  we  checked  the  t2  depen¬ 
dence  of  the  profiles  for  various  diffusion  times  (1.5,  6.0  and  24 
hrs  at  1000°C).  Regardless  of  the  diffusion  anneal  time,  all  profiles 
showed  well  defined  Regions  II  and  III  and  coincidence  with  one 
another.  We  conclude  that  Regions  II  and  III  are  bulk  diffusion 
controlled. 


Several  important  points  can  be  made.  These  are: 


(1)  Sulfur  and  tellurium  exhibit  virtually  identical  diffusion 
coefficients  in  both  modes  (Regions  II  and  III),  i.e.  both 
sulfur  and  tellurium  exhibit  the  same  slower  and  faster 
dif fusivities ,  Do  and  Df,  and  we  find  no  appreciable  size 
effect  [02"=1.32  A.  ,  S2“=  1.84  a,  Te2-=  2.11  A.  ]  or  mass  ef¬ 
fect  [0  18,  S  35,  Te  127]. 


(2) 


Both  modes  of  sulfur  and  tellurium  diffusion,  Ds  and  Df, 
exhibit  two  distinct  regions  in  the  temperature  dependence 
(See  Fig.  11).  At  the  higher  temperatures  (1050°C  ^T^  1250°C) 
both  modes  have  very  small  temperatures  (~  10  kcal/mole.) 

In  the  lower  temperature  range  (850°^Tsl050°C)  both  modes 
are  very  temperature  dependent.  The  Arrhenius  equations  for 
their  dif fusivities  in  cm2/sec  are 


850°^T^1050° ;  Ds  =  1.995  x  103  exp 
Df  =  6.84  x  10^  exp-[ 


[ -80.08  kcal/mol ’ 

-103 . 65^kca  1  /mo  1  - 
RT 


1050°^1125°C; 

Ds 


3.55  x  10 
3.09  x  10' 


exp  [“I?.3-2  l^l./moli 

ex d  r  ~f0. 25R£cal/mol i 
F  l  RT 


The  absolute  magnitudes  of  Ds  and  Df  for  S-35  and  Te-127  are 
almost  four  orders  of  magnitude  faster  than  S-35  in  NiO  reported  by 
Hwong  and  Wagner  (12).  This  result  is  in  general  agreement  with 
the  much  higher  concentration  of  oxygen  vacancies  in  CSZ  than  in  NiO. 
The  values  of  Ds  and  Df  are  less  than  the  diffusivitv  of  oxygen  in 
CSZ  but  greater  than  the  dif fus ivities  of  Ca  and  Zr. 
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Similar  experiments  were  carried  out  for  S-35  in  y ttria-stabilized 
zirconia,  both  polycrystalline  pellets  (9  m/o)  and  single  crystals 
(9.6  m/o).  In  these  studies  at  1000°,  1100°  and  1200°C  both  the 
solubility  of  sulfur  and  the  diffusion  coefficients  were  extremely 
small.  Single  and  polycrystalline  samples  yielded  the  same  values 
for  diffusivity  and  evidence  for  multimode  diffusion  was  not  obtained 
because  of  the  extremely  small  penetration  profiles.  Approximate  values 
were*.  .  ,  ~  , 

D(S-*YSZ )  s  8.68  x  10  cm  sec  (1200°C) 

2  1.16  x  10“14  cm2sec_1  (1100°C) 

-  3.78  x  10  ^  cm2sec  ^  (1000°C) 

In  related  experiments,  the  diffusivity  of  sulfur-35  was  deter¬ 
mined  in  CSZ  as  a  function  of  calcia  content.  Samples  of  9,  13,  15 
and  18  m/o  CaO  were  kindly  furnished  by  D.  Chatterji  and  K.  Lay 
of  the  General  Electric  Company. 

In  contrast  to  the  data  on  commercially  available  15  m/o  CaO- 
Zr02  from  Zircoa,  the  results  on  the  GE  material  depended  markedly 
on  the  time  of  the  pre-diffusion  anneal.  The  higher  the  pre¬ 
diffusion  anneal  temperature  and  longer  the  pre-diffusion  anneal  time, 
the  smaller  was  the  observed  sulfur  solubility  and  the  smaller  were 
the  diffusivities.  The  diffusion  profiles  so  obtained  did  not 
follow  the  expected  t'S  dependence.  Curiously,  the  electrical 
conductivity  showed  no  effects  of  "aging"  above  1100°C.  We  continue 
to  analyze  these  data. 
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Defect  Structure,  Ionic  Conductivity,  and 
Diffusion  in  Calcia-Stabilized  Zirconia 
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ABSTRACT 

A  quantitative  theoretical  model  which  accounts  lor  the  variations  ol  ionic 
conductivity  ol  calcia-stabilized  zirconia,  Ca*rZri-2iOj-£r,  with  oxygen  va¬ 
cancy  concentration,  x,  and  temperature,  T,  has  been  developed  based  on  the 
consideration  ol  effective  anion-site  coordination  ol  effectively  negatively 
charged  calcium  ion,  Ca2'.  The  theory  extends  the  conventional  concept  ol  dop¬ 
ant-defect  complexes  (or  associates)  and  demonstrates  that  ionic  conduction 
in  this  system  proceeds  by  a  multimode  mechanism  through  and  between  dil- 
ierent  degrees  ol  dopant-defect  associates.  The  primary  mechanism  is  oxygen 
transport  through  the  channel  of  onefold  effectively  Ca2'  coordinated  amon 
sublattice,  that  is,  oxygen  vacancy  Vo-  transport  inside  the  onefold  bonded 
(Ca2'-Vo2  )  defect  complex.  The  theory  subsequently  predicts  the  occurrence 
ol  a  rather  sharp  maximum  of  ionic  conductivity  at  oxygen  vacancy  concentra¬ 
tion  x  =  0.0825  and  reproduces  the  experimental  isothermal  log  <r  vs.  x  curves 
reported  by  various  investigators  quite  satisfactorily.  Accordingly,  numbers 
of  parameters  which  characterize  the  defect  structure,  the  ionic  conduction, 
and  the  diffusion  process  in  this  system  are  derived  from  the  curve- fitting  pro¬ 
cedure  of  the  theoretical  expression  for  ionic  conductivity  with  experimental 
ones. 


It  is  now  well  recognized  that  zirconia  stabilized  in 
the  cubic  fluorite  structure  by  substitution  of  aliovalent 
oxides,  such  as  GaO,  Y203,  and  rare  earth  oxides,  is  es¬ 
sentially  an  oxide  ion  conductor  with  its  high  ionic 
conductivity  due  to  high  oxygen  vacancy  concentration 
introduced  into  the  anion  sublattice  by  these  substitu¬ 
tions  (1).  Numerous  experimental  works  on  electrical 
conductivity  on  these  systems  have  shown  that  ionic 
conductivity  of  these  systems  exhibits  a  maximum 
value  around  the  minimum  dopant  level  required  to 
stabilize  the  cubic  fluorite  phase  at  constant  tempera¬ 
ture.  For  example,  in  CaO-stabilized  zirconia. 
CaizZrt-hOi-iz,  the  maximum  occurs  around  12-13 
mole  percent  (m/o)  CaO  content  (2-5)  (x  =  0  060-0.065) 
and  in  the  Y203-stabilized  zirconia.  Y<I2ri-<I02-2I, 
around  8-9  m/o  Y203  content  (2,4,6)  (x  =  0.037- 
0.041).  At  higher  concentrations  within  the  cubic  fluo¬ 
rite  phase  regions,  the  ionic  conductivity  decreases 
with  increasing  dopant  content.  The  activation  energy 
for  ionic  conduction  is  also  reported  to  increase  with 
increasing  dopant  content  i4,  5).  Similar  behavior  is 
also  reported  in  the  other  fluorite  oxide  systems,  i.e., 
Th02  (7)  and  Ce02  (8)  based  solid  electrolytes,  where, 

*  Electrochemical  Society  Active  Member. 

Key  won!*:  electrolyte,  conductance,  transport. 


since  pure  thoria  and  ceria  have  cubic  fluorite  struc¬ 
ture,  dopants  such  as  CaO  and  Y2Os  are  used  only  to 
introduce  the  extrinsic  oxygen  vacancies  into  the  anion 
sublattice.  This  decrease  in  ionic  conductivity  with  in¬ 
creasing  vacancy  concentration  has  been  qualitatively 
regarded  as  the  result  of  vacancy  ordering  (5)  va¬ 
cancy  clustering  (9),  or  dopant- vacancy  association 
(10)  in  the  higher  oxygen  vacancy  concentration  range. 

Recently,  several  attempts  to  explain  this  behavior 
have  been  undertaken:  O’Keeffe  ill)  and  Barker  and 
Knop  (12)  have  shown  that  when  vacancy-vacancy  re¬ 
pulsive  interaction  is  so  strong  that  the  numbers  of  the 
first  and  the  second  nearest  neighbor  vacancy-vacancy 
configurations  are  prohibited,  the  maximum  in  tome 
conductivity  shifts  to  the  lower  vacancy  concentration 
ranges  than  that  expected  from  the  usual  hopping 
model  in  which  v  ac  xl  1  —  x),  so  that  is  expected 
at  x  =  0.5.  Casselton  (6)  for  the  Y203-Zr02  system 
and  Hammou  (13)  for  the  YiC^-ThO;  system,  respec¬ 
tively,  proposed  models  to  explain  this  behavior  using 
the  extended  Lidiard  model  (14)  which  is  based  on 
the  presence  of  dopant-vacancy  complexes  and  Debye- 
Hiickel  type  interactions.  Hammou  also  tested  Barker’s 
model  mentioned  above.  None  of  these  models  resulted 
in  quantitative  agreement  with  experimental  data. 
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Most  recently,  Schmalzned  (13)  Interpreted  these  phe¬ 
nomena  under  the  assumptions  of  (t)  Debye-Huckel 
type  electrostatic  interaction  tor  low  detect  concen¬ 
tration  range  and  (ii)  ordering  o f  oppositely  charged 
delects  tor  high  detect  concentration  range,  neglecting 
the  formation  of  detect  associates.  Schm aimed  showed 
that  these  assumptions  lead  to  a  decrease  ot  the  cor¬ 
relation  factor  of  delects  with  increasing  delect  con¬ 
centration  and  also  to  an  increase  in  the  activation  en¬ 
ergy  for  delect  motion  in  a  qualitative  manner.  Judg¬ 
ing  from  the  preceding  brief  literature  survey,  there 
seems  to  exist  no  theoretical  model  which  accounts  for 
the  phenomena  in  these  systems  in  a  quantitative 
manner.  Accordingly,  the  aim  ot  the  present  paper  is  to 
propose  an  alternative  theoretical  model  to  serve 
this  purpose  by  choosing  calcia-stabilized  zirconia  as 
a  representative  system. 

In  the  next  section,  a  brief  survey  of  the  ionic  con¬ 
ductivity  data  of  calcia-stabilized  zirconia  reported  in 
literature  is  given,  and  in  the  subsequent  section  the 
theoretical  model  is  described. 

Ionic  Conductivity  Data  of  Calcia-Stabiiized  Zirconia, 
CauZn-trOj-tE,  System 

Figure  1  shows  the  isothermal  log  a  vs.  oxygen  va¬ 
cancy  concentration  x  plots  at  T  =  800*,  800‘,  1000*, 
1200*,  and  1400* C  reported  in  literature.  Included  are 
data  only  by  the  authors  who  studied  a  as  a  function 
of  dopant  content  in  a  systematic  manner. 

The  data  of  Strieker  and  Carlson  (2),  Tien  (3), 
Dixon  et  al.  (4),  and  Tien  and  Subbarao  (3)  all  agree 
that  the  stability  region  of  the  cubic  fluorite  phase  ex¬ 
tends  from  around  12-13  m/o  CaO  to  around  20-22  m/o 
CaO  at  temperatures  between  approximately  800*  and 
1400*  C.  Furthermore,  the  conductivity  maximum  exists 
inside  the  cubic  fluonte  phase  near  the  monoclinic- 
cubic  solid  solution  phase  boundary,  t.e.,  at  12-13  m/o 
CaO  content  tx  =  0.080-0.085).  In  contrast.  Carter 
and  Roth  <  18)  report  that  the  stability  region  extends 
from  about  10  to  19  m/o  CaO  at  1400'C  and  their  data 
seem  to  show  the  tendency  that  the  conductivity 
maximum  shifts  from  around  13  m/o  CaO  at  1000*C  to 
15  m/o  CaO  at  1400*C. 

As  is  apparent  from  Fig.  1,  the  ionic  conductivity  of 
calcia-stabilized  zirconia  initially  increases  with  oxy¬ 
gen  vacancy  concentration  x.  and  reaches  a  maximum 
around  12-13  m/o  CaO  content  .  -  =  0.080-0.085).  The 
conductivity  then  decreases,  initially  rather  sharply 
and  then  more  slowly,  with  oxygen  vacancy  concentra¬ 
tion  x  within  the  cubic  fluonte  phase  region.  The  rate 
of  decrease  of  ionic  conductivity  with  x  increases  with 
decreasing  temperature  and  the  conductivity  decreases 
by  a  factor  of  —  1/20  at  800*  C  when  r  increases  from 
0.065  to  0.100.  This  means  that  the  activation  energy  for 
ionic  conduction  increases  with  x.  Typical  values  are 
28.2  Iccal/mole  for  13  m/o  CaO  and  31.1  kcal/mole  for 
20  m/o  CaO  (5). 

Another  interesting  feature  is  the  order-disorder 
transition.  Subbarao  et  al.  1 5,  171  observed  this  phe¬ 
nomenon  below  1000'C  in  samples  with  higher  CaO 
content  (18  and  20  m/o  CaO).  Some  of  these  data, 
designated  as  "ordered,'*  are  also  included  in  Fig.  1. 
"Ordered”  samples  at  these  compositions  have  a  lower 
but  finite  conductivity  which  is  converted  to  that  of 
the  'disordered"  phase  after  heat-treatments  at  1400*C. 
Carter  and  Roth  (18)  also  examined  this  order-dis¬ 
order  phenomenon  rather  extensively.  Their  data  in¬ 
dicate  that  below  ilOO’C,  even  in  samples  with  low 
CaO  content  ( 13.2  and  14.2  m/o  CaO),  in  some  cases 
the  resistivity  changes  by  a  factor  of  over  10.  show¬ 
ing  no  tendency  to  saturate  with  time.  The  reason  for 
this  discrepancy  between  Subbarao  et  al.  (5.  17)  and 
Carter  and  Roth  1 16)  is  not  clear  at  present,  although 
it  seems  reasonable  to  assume  from  these  observations 
that  in  samples  with  higher  CaO  content  and  at  lower 
temperatures  calcium  ions  Ca2'  and  oxygen  vacancies 
Vo3  would  order  in  some  fashion  which  causes  the 


!V02  l  =  x 

Fig.  I.  Iiotbormal  log  a  rt.  oxygon  vacancy  concentration  < 
plats  at  600‘,  800*.  1000’,  1200°,  and  1400°C 

Reference 


Experimental  points  number 

C  Strickler  and  Carlson  (2) 

A  Tien  (3) 

X  Dixon,  LaGrange,  Merten,  Miller  and  Porter  (4) 

V  Tran  and  Subbarao  (5) 

^  Carter  and  Rath  (16) 


—  Theoretical  curves  according  to  Eq.  [13]  at  0.0625 
=si  x  0.125  and  Eq.  [13]  at  0  =£  r  0.0625,  T  = 
1000°C 

[  The  phase  boundary  of  the  cubic  fluorite  phase  ac¬ 
cording  to  the  above  authors. 

VAovrer  values  indicated  by  arrow,  a  at  ordered  state 
according  to  Ref.  (5)  and  (17). 

substantial  decrease  in  ionic  conductivity  after  pro¬ 
longed  heat- treatments. 

Theoretical 

The  main  aim  of  the  theoretical  model  is  to  repro¬ 
duce  the  experimental  log  a  vs.  x  curves  shown  in  Fig. 
1  over  the  entire  temperature  ranges  in  a  quantitative 
manner  consistent  with  the  physical  reality  of  the  sys¬ 
tem.  This  includes  the  increase  in  activation  energy  for 
ionic  conduction  with  x. 

As  far  as  the  authors  are  aware,  Kroger  (10)  was  the 
first  to  suggest  that  ionic  conduction  and  oxygen  dif¬ 
fusion  in  calcia-stabilized  zirconia  (CSZ)  probably  oc¬ 
curs  via  migration  of  oxygen  vacancies  Vo2,  in 
(CazrV'o)-/.  Kroger  described  this  association  process 
by  a  usual  defect  equilibrium  equation 

mCaZr2’  4-  mV03  ^  (CaZrV0)m*  [1] 

Applying  the  mass  action  law  to  the  above  equation 
and  assuming  electrostatic  binding  energy  of  the  defect 
complex.  Kroger  concluded  that  most  of  the  calcium 
(and  therefore  Vo2  )  is  present  as  neutral  calcium  ion- 
oxygen  vacancy  pairs  or  clusters.  This  description  cor¬ 
responds  in  general  form  with  the  theory  of  Lidiard 
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tor  aliovalent  impurity  diffusion  in  alkali-halide  sys¬ 
tems  (14).  But  when  the  impurity  ( doper  t)  concentra¬ 
tion  and  therefore  also  the  delect  concentration  be¬ 
comes  as  high  as  in  CSZ,  these  treatments  do  not  ap¬ 
pear  to  be  valid.  First,  let  us  assume  that  oxygen 
and  therefore  Vo3'  diffusion  proceeds  vis  delect  com¬ 
plexes  (CazrVo)1  (m  s  1  in  Eq.  [1]).  In  the  fluorite 
structure,  since  every  cation  has  8-flrst,  24-second,  24- 
third  .  .  .  nearest  neighbor  anion  sites,  all  the  anion 
sites  will  be  just  effectively  onefold  coordinated  by 
dopant  cations  Ca3’  when  the  dopant  concentration 
reaches  2x  =  UZ'  tor  Cai*Zri-iiOt-tn  i-e.,  in  this  case 
Z’  -  2x  =  2,  where  Z'  is  the  effective  anion  site  coordi¬ 
nation  number  ol  Ca3’  in  this  system.  This  means  that 
at  the  composition  r  =  1/Z’,  every  V03-  jump  occurs 
inside  (Ca3’-V03  )  onefold  effectively  associated  delect 
complexes.  Therefore,  conductivity  should  be  a  maxi¬ 
mum.  Beyond  this  composition  the  average  effective 
Ca3’  coordination  number  to  the  anion  sublattice  be¬ 
comes  greater  than  1,  i.e.,  the  number  ol  density  ol 
2-fold  Ca3’  coordinated  anion  sites  becomes  an  ap¬ 
preciable  fraction  ol  the  total  available  jump  sites; 
so  the  ionic  conductivity  should  decrease,  because  the 
effective  path  for  Vo3'  migration  through  (Ca2’-Vo2'): 
1-fold  effectively  Ca3’  associated  defect  complexes 
decreases.  This  is  indeed  the  case  lor  CSZ  if  we  assume 
an  effective  anion  site  coordination  number  ol  Ca3’:  Z’ 
=  16:  the  experimental  observed  ionic  conductivity 
maximum  occurs  around  x  =  0.060-0.065,  whereas  the 
theoretically  derived  conductivity  maximum  is  ex¬ 
pected  at  x  =  1/Z’  =  1/16  =  0.0625.  If  only  the  flrst 
nearest  neighbor  Ca^-Vo3'  configuration  is  taken  into 
account,  Z’  becomes  equal  to  8,  and  the  conductivity 
maximum  is  expected  at  x  =  1/8  =  0.125,  which  is 
almost  twice  the  observed  oxygen  vacancy  concen¬ 
tration  at  the  maximum.  This  means  that  the  second 
and  the  higher  order  nearest  neighbor  Ca3’-Vo3'  con¬ 
figurations  also  contribute  to  the  effective  anion  site 
coordination  number  ol  Ca3’,  Z’,  with  a  contribution 
that  is  equal  to  that  from  the  first  nearest  neighbor 
0a3’-Vo3  configuration  (Z*  =  8).  This  fact  seems  to  be 
quite  reasonable  in  the  light  of  the  rather  long-range 
character  ol  electrostatic  interaction  between  charged 
delects  Ca3’-V03'  (18).  We  might  regard  the  first  near¬ 
est  neighbor  Ca3’  and  Vo3-  configuration  as  the  ground 
state  ol  the  (Ca3'-V03)  defect  complex,  and  the  second 
(and  the  higher  order)  nearest  neighbor  Ca3'-Vd3' 
configuration  as  the  first  excited  state  ( and  the  higher 
order  excited  states)  ol  the  (Ca3'-V03  )  defect  complex. 

But  actually  the  greater  the  distance  between  Ca3' 
and  V02  .  the  more  the  electrostatic  interaction  be¬ 
tween  them  would  be  increasingly  screened  and  weak¬ 
ened.  Consequently  we  include  here  only  the  first  near¬ 
est  neighbor  and  the  second  nearest  neighbor  Ca3'- 
V03  configurations,  assuming  further  that  Ca3’  coordi¬ 
nation  number  to  each  24-second  nearest  neighbor 
anion  sites  is  —  1/3.  At  this  level  of  approximation, 
the  total  effective  coordination  number  of  Ca3’  to  anion 
sites  Z’  becomes  equal  to  (8  x  1)  +  (24  x  1/3)  =  16, 
which  yields  the  conductivity  maximum  at  x  =  0.0625 
in  accordance  with  experimental  observations. 

When  the  CaO  content  becomes  25  m/o,  i.e..  x  = 
0.125,  the  effective  anion  site  coordination  number  of 
Ca3'  becomes  2.  i.e.,  Z’  ■  2x/2  =  2.  This  situation  is 
shown  in  Fig.  2a.  Figure  2a  shows  the  unit  cell  of 
Cao,2jZro,7jOi.7jVoo.«  (x  =  0.125).  The  unit  cell  of  the 
fluorite  structure  contains  4  cation;  and  8  anions. 
So  at  the  composition  2x  =  1/4,  one  out  of  four  cations 
is  Ca3’.  In  Fig.  3  these  Ca3'  are  situated  at  the  eight 
corner  cation  sites  of  the  unit  cell  The  remaining  three 
cations  are  Zr1.  which  are  situated  at  the  face  center 
position  of  six  face  planes.  One  out  of  8  anions  in  the 
unit  cell  is  V02  (x  =  1/8).  As  is  apparent  from  this 
figure,  i/4  of  the  four  (=1)  first  nearest  neigh¬ 
bor  cations  around  V02  is  Ca3’,  which  yields  onefold 
Ca3'  coordination  to  V02\  and  1/4  of  the  twelve  (=  3) 
second  nearest  neighbor  cations  around  V02-  are  Ca3’, 


Fig.  2s.  Unit  call  of  ths  fluorite  lattice  of  (CaO)o  2s(ZrOs)o  rs 
(2*  =  1/4,  x  =  0.125). 


Fig.  2b.  Ca2'-Vo3  configuration  in  (CoO)o  i2s(ZrOt)o  875(2*  = 
1/8,  »  =  0.0625).  •  Ca3',  0  Zr1,  G  O*.  ■  Vo3.  —  4  first  near¬ 
est  neighbor  cations  around  V0-  ,  — .  3  lout  of  12)  second  nearest 
neighbor  Co2'  around  Vo3  • 

whch  yields  another  onefold  coordination  to  V03'.  So, 
at  this  composition  every  V03'  (and  anion  site)  is  on  an 
average  twofold  Ca3’  coordinated  and  every  V03  jump 
occurs  inside  effectively  twofold  associated  (Ca2’-Vo3  - 
Ca3’)  defect  complexes.  Similarly,  Fig.  2b  shows  the 
Ca3’-Vo2  configuration  in  Cao.12sZro875Ot.875Voo.12s2' 
(x  =  0.0625)  where  every  two  unit  cells  contain  one 
Ca2’  (2x  =  1/8)  and  one  Vo2  (x  =  1/16).  In  Fig.  2b. 
Ca3’  are  situated  at  face  centered  positions  of  the 
top  plane  or  bottom  plane  in  an  alternative  array. 
Every  V0-  (and  anion  site)  is  just  onefold  coordinated 
by  Ca- ,  either  by  one  first  nearest  neighbor  Ca3’  ( =  1) 
or  by  three  second  nearest  neighbor  Ca3’  ( 1/3  x  3 
=  1).  So  at  this  composition  every  Vo3  jump  occurs 
inside  effectively  onefold  associated  (Ca3’-Vo3)  defect 
complexes. 

Note  that  the  stabilized  cubic  fluorite  phase  appears 
just  when  onefold  associated  <Ca3’-V03  )  defect  com¬ 
plexes  extend  throughout  the  anion  sublattice,  giving 
rise  to  a  conductivity  maximum  (x  =  0.0625) ;  and  dis¬ 
appears  with  the  disappearance  of  th,s  defect  complex 
<x  =  0.125).  In  these  defect  concentration  ranges,  con¬ 
ventional  defect  equilibria  treatments  such  as  Eq.  [1] 
are  not  applicable,  for  there  exist  no  more  free  dissoci¬ 
ated  anion  sites  on  which  Vo-’  can  reside. 

Description  of  the  Theoreficol  Model 

Let  us  designate  at  a  given  composition  x,  the  rela¬ 
tively  free  state,  and  the  relatively  bound  state  of  the 
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anion  sublattice  tor  Vo*'  aa  F  and  B,  respectively.  In 
the  case  of  CSZ  within  the  vacancy  concentration  range 
0.0629  *  x  *  0.129  over  which  the  cubic  fluorite  phaae 
extendi,  we  can  that  the  F-state  and  the 

B-atate  correspond  to  the  1-fold  Ca*'  coordinated 
and  2-fold  Ca**  coordinated  anion  sublattice,  respec¬ 
tively.  Thii  follows  because  the  Ca*'-average  eoordi- 

Z*  •  2* 

nation  number  to  the  anion  sublattice  is  1  ^  — - — 

2 

—  isx  ^6  2.  Tor  each  concentration  of  E-state  and 
B-state,  the  site  tractions  Nr  and  NB  will  be  given  by 
the  following  equation 

1  x  N»  +  2  X  N,  =  Z’  ■  3x/2  =  16*.  Nr  +  NB  =  L 

Nr  *  2  -  18*.  Na  =  16*  -  1 

0.0629  **«S  0.129  [2] 

Equation  [2]  neglects  the  statistical  fluctuations  of  the 
system.  The  distribution  of  oxygen  vacancies  Vo*'  and 
oxide  ions  O*  between  these  states  is  assumed  to  be 
given  by  the  following  quasi-chemical  equilibria 

Vo*-  (B)  +  0*<F)  v*  Vo*-  (F)  +  0»(B)  [3] 

Applying  the  mass  action  law  to  Eq.  [3]  yields 

K  =  e*p(-^l)=e*p(i£l)  .exp(-^l) 

_  [Vo*(F)3lO«(B)] 

"  [Vo*<B)][0*(F)] 

where  K,  aG‘,  aH‘,  and  denote  the  equilibrium 
constant,  standard  Gibbs  free  energy  change,  standard 
enthalpy  change,  and  standard  entropy  change  for  Eq. 
[3],  respectively,  and  [  ]  designates  the  concentration 
of  each  chemical  species  in  site  fractions. 

From  the  site  and  mass  conservation,  the  following 
equations  hold 

[Vo*- (total)  J  =  [Vo*(F)]  +  [Vo*(B)J  =* 

Nr  +  N,  =  1 

[Vo*(F)I  +  tO*(F)  ]  =Nr 

[Vo*(B)]  +  [0*(B)]  =Na  [9] 

Rearranging  Eq.  [9]  yields 

[Vo*  (B)J  =  *-  [Vo*- (F)] 

[0*(F)1  =Nr-  [Vo*  (F)] 

[0*(B)]  =  1  —  Nr  —  *  +  [V0(F)]  [6] 

Inserting  Eq.  [6]  into  [4],  the  following  quadratic 
equation  for  [V0*-  (F)  ]  is  obtained 

(1  -K)[Vo*(F)]*  +  {K(*  +  Nr) 

+  a-Nr-*)}[Vo*  (F)J  -K  Nr-x  =  0  [7] 

Solving  Eq.  [7]  for  [Vo*  (F)]  (>  0)  yields 


f  / 

c  TO"1  :  / 

b  :  / 


IV0J1»  x 

Fig.  1.  Schematic  variatiox  of  »  with  oxygen  racancy  concmtra- 
Hon  x  whan  Vo*'  transport  only  through  Co2’- 1 -fold  coordinated 
onion  inblaftki  it  attuned.  Number?  betide  each  caret  at  0.0625 
did 0.125  doi«Mta  K  ralssi  of  Eq.  [16], 

the  channels  of  the  1-fold  Ca*'  coordinated  and  2-fold 
Ca*'  coordinated  anion  sublattice,  respectively,  and 
<r  and  <t  designate  those  due  to  V0*-  transport 

Ih2  2dl 

from  1-  to  2- (or  reverse) fold  Ca*' -coordinated  anion 
sublattice. 

Each  partial  ionic  conductivity  will  be  given  by  the 
following  expression  utilizing  the  usual  form  for  hop¬ 
ping  ionic  conduction 


N  Z-e 


=  [V0*(i>] 


[Vo*  (3)] 


[V0*  (i)](0*(j)] 


N,  +  Nj 

where  ««>  denotes  the  mobility  of  Vo*-  for  the  i-j 
jump  process.  Each  partial  ionic  conductivity  o't-j  is  nor¬ 
malized  by  the  numerical  factor,  (N  ■  Z  •  e)~l  from 
the  actual  ionic  conductivity  #vj.  In  this  expression  N  is 
numbers  of  an ion  sites  per  unit  volume  (cm3),  Z  is 
valence  number  of  Vo*  ,  and  e  is  the  electronic  charge. 
From  dynamic  equilibria  between  1-  and  2-states 


a  ss  a 

l-»i  3-.1 


Combination  of  Eq.  [4],  [10],  and  [11]  yields 


'*'*-«>  /  ag*  \ 

— =exp(-Rr) 


[Vo*  (F)]  = 


-{K(*  +  Nr)  +  (l-Nr-*)}  +  V(K(i 

2(1 


Then,  from  Eq.  [8]  and  [6],  all  the  concentrations  of 
each  chemical  species  can  be  obtained  given  K  and  r. 

The  total  ionic  conductivity  of  the  system 
will  be  given  by  the  sums  of  partial  ionic  conductivity 
»H  (i,  j  =  F,  B) 


..  2  n*.  X 


For  simplicity,  the  indexes  F  and  B  were  designated  by 
1  and  2  respectively,  where  »i.i  and  n-i  designate  the 
partial  ionic  conductivity  due  to  Vo*-  transport  through 


+  Nr)  +  (1  -  N,  -  *)}*  +  4*K(1  -  K)Nr 

-K)  - -  181 

Equation  [12]  is  the  well-ltnown  relationship  for  chem¬ 
ical  equilibria  resulting  from  opposing  elementary 
reaction  processes. 

Inserting  Eq.  [10]  into  Eq.  [9],  the  Anal  expression 
for  the  normalized  total  ionic  conductivity  of  the  sys¬ 
tem  is  obtained 


9  norm«Ux«l  =  * 


—  =  <r'i.i  +  /J4+  2<r'^j 


=  [Vo*-(l)][0*U)Wn  +  [Vo*  (2)][0*(2)]M,M, 


+  2[Vq*  (2)] [O1  (!)]/» 


(N,  +  N,=  1)  [13] 
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where  [V02(i)]  and  CO*(i)]  (i  =  1,  2) are  given  by 
Eq.  Lfl]  and  [8]  tor  a  given  K  and  x. 

Bosk  Character  of  Theoretical  Model 
We  will  first  examine  the  theoretical  model  in  more 
qualitative  but  illustrative  manner.  For  this  purpose, 
tar  mathematical  simplicity,  we  will  adopt  the  follow¬ 
ing  approximation  in  this  section:  [Vo2'(i)J  << 
[OMili  m  Nt  li  =  F,  B),  (which  corresponds  to  ne¬ 
glecting  the  (1  —  x)  term  in  the  conventional  a  a 
x(l  —  x)  form  by  assuming  (1  —  x)  —  1).  At  this 
level  o t  approximation,  Eq.  13]  can  be  easily  solved  to 
give 

x  •  Nr  •  exp  (-uGVRT) 

[V0,-(F)]=  IVa^r-expl-nG-B-r 

x  •  Nr  •  K 
“  N9  +  Nr  ■  K 


a’ noro*U»d<iauJ  =  »’l-l  (x  s  0.0825)  —  0. 0825m m>  [17] 

However,  in  this  approximation,  according  to  Eq.  [16] 
ionic  conductivity  becomes  equal  to  zero  at  x  =  0.125 
where  the  Ca2'-l-fold  coordinated  anion  sublattice 
disappears.  At  this  composition,  only  the  Ca2’-2-fold 
coordinated  anion  subiattice  exists  (see  Fig.  2a),  so 
from  Eq.  [14]  we  obtain 

a  Qonujiild  (St  —  0.125)  z  a  2-3  (X  —  0.125) 

=  0.125 -wsa)  [18] 

Equations  [17]  and  [18]  indicate  that  from  the  con¬ 
ductivity  data  at  x  =  0.0625  and  0.125,  we  can  obtain 
the  mobility  mhi  and  mm)  inside  the  separate  chan¬ 
nels.  Using  these  values  of  mi  d,  ms-2>.  we  can  draw  the 
theoretical  log  a  vs.  x  curves  over  the  whole  range  of 
0.0625  s  x  s  0.125,  choosing  appropriate  values  of  K 


[TV  (B)]  =  x-  [Vo»-(F>] 
*'N, 


x-N  a 


NB  +  Nr  ■  exp  (— hGVKT)  NB  +  Nr  ■  K 


[8-1 


Using  Eq.  [2]  for  Nr  and  Nb,  and  also  Eq.  [13]  for 
- - — i  the  following  expression  for  v'aonnaiiaad  is  ob¬ 
tained,  valid  for  0.0625  ^14  0.125 


#tocai 


N  Z-t  (18x  -  1)  +  (2  -  16x)  ■  K 
X  {(2  -  ldx)2  •  K  •  mu)  +  2(2  -  18x)  (16x  -  1) 

•  +  (18l  —  1)2M14>}  [14] 


Numerical  Evaluation  of  the  Theoretical  M'del 
In  order  to  obtain  mobilities  mm>  and  ms-*>.  ap¬ 
proximate  values  of  vault  ous)  —  ai-t  (x  =  0.0025),  and 
VtOUltX  —  0.  125)  =  =  0.125)  were  read  from 

experimental  log  a  vs.  x  plots  shown  in  Fig.  1.  For  x  = 
0.125  which  is  well  beyond  the  cubic  fluorite  phase  so 
that  experimental  points  for  awui  are  not  available, 
extrapolations  from  the  data  points  around  x  ~  0.100 
composition  were  made.  Moreover,  the  values  of  >au 
at  x  =  0.625  and  a wtu  at  x  =  0.125,  were  also  adjusted 
to  get  the  best  fit  between  theoretical  and  experi¬ 
mental  curves  in  the  subsequent  calculations.  At  a 
given  pair  of  such  »’s,  mm>  and  mm>  were  calculated 
irom  isq.  [17' J  and  [18'] 


Following  the  same  procedure,  for  the  composition 
range  0.0625 

N,  s  No  =  (1  -  16x) :  (CaJ'-0-fold  state) 


a  max  =  ff'i-i(x  =  0.0625)  2  X(1  - 


■  x)  •  m-i> 

15 

= — •  «i-n  [in 


N,  =  N,  a  16x:  (Ca2'-l-fold  state) 

and 

x 

'mdM-  ldx  +  (1  —  16x)  ■  K’  *  ^ 


18x)2  ■  K’ 


a  coui  (x  —  0.125)  ns  a  ss(x  —  0,125) 


=  X(1  —  X)  ■  MM)  2 


T_ 

8* 


'  WM) 
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'  M<M>  +  2  •  (1  —  ldx)  •  18x  +  (16X)2  •  Mi'll} 

[15] 

w^ere  K'  =  exp  (— dG'VRT)  and  uG“  designates  the 
Gibos  free  energy  difference  between  Ca2'-0-fold  and 
1-fold  states. 

When  aG-  =  aG'  =  0  (K’  =  K  =  1),  all  the  sites 
become  energetically  equivalent,  so  all  the  mi.j>  be¬ 
come  equal  and  .h._ <  «  x  from  Eq.  [14]  and  [15], 

Another  interesting  limiting  case  is  the  one  in  which 
ionic  conduction  occurs  exclusively  inside  the  chan¬ 
nel  of  the  Ca2'- 1-fold  coordinated  anion  sublattice, 

(i.e.,  K’  <<  1,  MM)  >>  MM).  K  <  1,  MM)  >> 

Mt-ai,  it  ).  In  this  case  from  Eq.  [14]  and  [15] 

<j*»n 

^'normalis'd  =  *  '  18X  '  MM),  0  ^  X  ^  0.0625 
x(2  —  16x)2  MM)  •  K 
(18x  -  1)  +  (2  -  16x)  •  K 

0.0825  ^x^0.125  [16] 

Figure  3  shows  /nonuibu/wM)  vs.  x  curves  using 
K  =  exp  (— sG’/RT )  as  parameter.  Included  also  is 
the  usual  curve,  #  «  x.  It  is  apparent  from  this  figure 
that  if  ionic  conduction  proceeds  only  via  Ca2’- 1-fold 
associated  (Ca2'-V02  )  defect  complexes,  there  appears 
a  rather  sharp  conductivity  maximum  at  x  =  0.0625, 
and  then  a  decreases  more  rapidly,  the  smaller  the 
value  of  K. 

These  features  are  qualitatively  in  fair  agreement 
with  experimental  observations  as  shown  in  Fig.  1. 

From  Eq.  [16]  the  maximum  conductivity  at 
x  a  0.0625  is  given  by 


Using  these  values  of  mm>  and  mm>,  and  assuming  ap¬ 
proximate  values  of  K  and  mm>,  the  theoretical  log 
a  vs.  x  curves  were  constructed  at  each  temperature 
in  intervals  of  ax  =  0.005  for  0.0625  as  x  &  0.125  until 
the  best  fit  to  the  experimental  curves  was  obtained. 

At  T  =  600’  and  800 ‘C,  the  best  fits  are  already  ob¬ 
tained  without  introducing  »  terms,  i.e.,  assuming 

a  —  0  ■  •  <r  —  0.  At  higher  temperatures  T  = 
2bi 

1000",  1200",  and  1400"C.  the  relative  contribution 
from  the  partial  ionic  conductivity  a  becomes  ap¬ 
preciable,  and  to  obtain  the  best  fit  these  contributions 
cannot  be  neglected. 

Results  and  Discussion 

The  best  fit  theoretical  curves  obtained  using  the 
rigorous  formula  for  awui  (Eq.  [13])  and  for  concen¬ 
trations  of  the  respective  chemical  species  as  given  by 
Eq.  [6]  and  [8]  are  shown  in  Fig.  1  as  solid  lines  at 
each  temperature  within  the  composition  range  0.0625 
^  x  ^  0.125.  Numbers  of  characteristic  parameters  used 
to  construct  these  curves  are  tabulated  in  Tables  I  and 
II.  .Also  included  is  the  theoretical  curve  for  aloui  for 
0  ^  x  0.0625  at  100(TC  when  it  is  assumed  that  only 
#m  contributes  to  the  ionic  conductivity:  i.e. 

a M  _ 

a  normalis'd  —  -  —  X(10X  —  X)  *  MM) 

♦V  •  Z  •  e 

ssx-lSx-mi.it,  0  ^  x  —  0.0625  [13’] 

Agreement  between  the  theoretical  curves  and  the  ex¬ 
perimental  ones  is  quite  satisfactory  over  the  entire 


i 


,c 
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Tabid  I.  K,  I  mu  U  =  0.0623),  and  a  (*  =  0.125)  mod  la 
com  true  I  Hm  Hi  tor*  heal  cutyw  in  Fig.  1 


T  CC) 

d«u  ™  6i-i 

(X  -  0.0825) 

lU  •  cm)-4 

<Td-l  m  j 

(x  -  0.125) 
lU  cm)-4 

K 

600 

7.341  x  10- 

5.203  x  10-* 

0.080 

800 

1.165  x  10- 

4.0  x  10“* 

0.12 

1000 

7.015  x  10-* 

4.3  x  10- 

0.20 

1200 

2.70  x  10-‘ 

4.2  <  10- 

0.26 

1400 

6.40  x  10-4 

2.0  v  10-  ‘ 

0.37 

Table  II.  Mobility  at  orygan  vacancy  Vo2  mvo2  (■  —  j)  (i  =  1.2) 
tbrmgb  and  between  Co2'- Mold  and  Ca2'-2-fold  coordinated 
mm  leblettke 

T 

fCi 

ItU-U 

[Cm*/ 
sec  V) 

Mte-» 

(cm*/ 

sec  V) 

Mi*-.. 

icm*. 

sec  V) 

M.,_g. 

( cm *< 
sec  V) 

600 

6.684  x  10-? 

2.507  x  10- 

8U0 

1.048  x  10- 

1.827  x  10-- 

_ 

1000 

6406  x  10-- 

2.168  x  10- 

6.300  x  10- 

343  a  10- ■ 

1200 

2.428  x  10  * 

2.024  <  10- 

1.500  x  10- 

5437  x  10-- 

1400 

5.758  x  10-* 

9.640  x  10- 

7.300  x  10- 

2.027  x  10-* 

temperature  range:  600*C  ^  T  ^  1400'C  and  the  sta¬ 
bility  range  of  the  cubic  fluorite  phase:  0.0623  ^  x 
=s  0.10. 

Figure  4  shows  Arrhenius  plots  of  #u(i)  (i  =  1,  2; 
x  =  0.0625,  0.125)  and  <tm(x)  •  T  (i  =  1,  2;  x  =  0.0625, 
0.125).  From  this  figure,  we  can  obtain  the  Arrhenius 
equations  for  <n-i  and  «*u  -  T,  respectively 

=  0.0625)  =  <UI  =  1.138  -  1»> 

(  25.022  kcal/znole  \ 

exp  l - - )  (a '  cm,"‘ 


U00T  1200  ‘C  1000X  800 ‘C  600X 


Fig.  4.  lag  »i-,(»)  and  lag  ■  T  («  =  0.0625.  0.125,  i  = 

1.2)  ►».  1/r  plait. 


<rj. ttx  =  0.125)  s  2.317  •  10* 


da  i)  =  1.19  •  exp 


(- 


25.022  kcal/mole 
RT 


(  38.833  kcal/mole  \ 

■  exp  ^ - — - J  (a  ■  cm)  [19] 

and 

*i-i  •  ix  -  0.0625)  ■  T  =  ■  T  =  4.57  •  10* 

/  27.474  kcal/mole  \ 

eXPl - RT - /  (acm)"‘-sK 

<7i., (x  =  0.125)  •  T  =  7.18  •  107 


cm1  •  sec-1  •  V-* 


*7(2-2)  =  1.16  ■  10*  ■  exp 


( 


38.835  kcal/mole 
RT 


) 


Da-i )  =  1.77  • 


10” *  •  exp 


( 


cm2  •  sec-1  •  V-*  [23] 

27.474  kcal/mole  ^ 

RT  ) 


exp 


40.936  kcal/mole 
— 


) 


(a  •  cm)-1  •  °K  [20] 


To  obtain  the  actual  value  of  ionic  conductivity  the 
number  of  charge  carriers  per  cubic  centimeter  is 
necessary.  To  this  end.  the  average  lattice  constant: 
ao  =  5.131A  (1)  for  the  nominal  composition  Cao.is 
ZrosiOi  as  was  adopted,  neglecting  the  relatively  small 
variation  of  a0  with  composition.  Then 

«ulx)  ii  =  1.  2;  i  =  0.0625,  0.125) 

=  .V  •  Z  e  •  xll  -  x)  •  mu) 


=  N  ■  xil  -  x) 


&  ■  «»  •  D(U) 
kT 


[21] 


where  .V  is  the  number  of  anion  sites  per  cm3  and  given 
<V  =  5.922  x  1022/cmJ  using  the  above  a«  value.  In  the 
latter  part  of  Eq.  [211  the  Nernst-Einstein  relation  was 

used 


M(M)  = 


Z  ■  e  •  P(H) 

kT 


[22] 


Combining  Eq.  [19],  [20],  and  [21]  leads  to  the  A r* 
rhenius  equations  for  »#u>  and  D*u> 


cm3  •  sec-1 

_  (  40.936  kcal/mole  \ 

V  RT  J 

cm2  ■  sec-*  [24] 

Values  of  „<i.n  and  D(|.p  are  tabulated  in  Tables  II 
and  III,  respectively,  together  with  m< i-j>  (i  -»  j  =  1,  2) 
and  which  were  also  used  to  derive  the  theo¬ 

retical  curves.  To  obtain  u  from  <*  Eq.  [12] 

(2-»l)  L  ‘ 

and  K  values  in  Table  I  were  used. 

Figure  5  shows  Arrhenius  plots  of  the  various  dif¬ 
fusion  coefficients  of  the  oxygen  vacancy  Du-p.  From 
this  figure  it  is  apparent  that  D<i.»  >>  D^.,)  ~D 

Tabic  III.  Diffusion  coefficients  of  oxygen  vacancy  Vo*-  Dm-j) 
through  and  between  Co2'* Mold  and  Ca2’-2*fold  coordinated 
anion  sublattice 


T  Dm-,.  Du-*.  ^  -•* 

(*C)  <cm-sec>  icm-/sec)  cm-  sec*  i cm-  sect 


600 

2.514 

Y 

10- 

9.429 

X 

10-‘ 

_ 

800 

4.844 

X 

10- 

8.910 

X 

10- 

— 

_ 

1000 

3.459 

X 

10- 

1.189 

X 

10- 

3.565  x 

10- 

1.783  x 

10-' 

1200 

1.541 

X 

10- 

1.284 

Y 

10- 

9  521  x 

io -• 

3.400  x 

10-' 

1400 

4.149 

< 

10- 

b.946 

< 

10- 

5.407  < 

10-' 

1.461  * 

10- 

log  O’  i  -  i  I  ohm  *  cm1! 
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>  D  as  expected.  This  mean*  that  diffusion  of  Vo2- 

is-.ii 

through  the  relatively  free  state  (1-1)  proceeds  much 
more  easily  with  a  lower  activation  energy  and  with  a 
higher  absolute  value  than  that  through  the  relatively 
more  bound  state  (2-2). 

An  interesting  feature  of  D  and  D  is  the  fact 
(s*i>  u— s> 

that  there  exists  a  rather  sharp  kink  in  the  Arrhenius 
plots  for  O  and  D  around  and  below  1000*C. 

(S-l)  (1—3) 

Indeed,  as  was  mentioned  in  the  previous  section,  the 
theoretical  curves  at  T  =  600°  and  800"C  were  con¬ 
structed  without  any  contribution  from  a  ,  and  yet 

shows  good  agreement  with  the  experimental  ones.  At 
1000’C,  their  contributions  are  still  small  but  at  higher 
temperatures  T  =  1200°  and  1400°  C,  they  become  con¬ 
siderable.  It  seems  premature  to  draw  any  definite 
conclusion  concerning  the  origin  of  this  phenomenon. 
However,  it  is  interesting  to  note  that  even  in  “dis¬ 
ordered”  samples,  it  is  often  reported  that  there  seems 
to  exist  a  slight  hump  in  the  Arrhenius  plots  of  log 
a  vs.  1/T  above  1000'C,  which  becomes  more  apparent 
in  the  samples  with  higher  CaO  content  (1,  5).  Figure 
6  shows  the  log  K  vs.  1/T  plot  for  the  chemical  equi¬ 
libria  given  by  Eq.  [3]  and  tabulated  in  Table  I.  A 
good  straight  line  relationship  holds  for  600°  C  T 
ss  1400°C.  From  this  figure  we  obtain  for  the  equilib¬ 
rium  constant  K 


/  AS'  \ 

/  AH'  \  I 

AG'  \ 

K  =  exp(_) 

.  ~  RT  ) 

8.606  kcal/mole  ^ 

L4’] 

—  2.7  ■  exp  (  - 

RT  ) 

Accordingly 

aS°  =  0.993  eu  = 

1.974  cal/mole  °K~R 

iff’  =  6.606  kcal/mole 

iG*  =  iH”  -  TaS’  =  8.608  x  103  -  1.974  •  T  cal/mole 

[23] 

Figure  7  shows  the  variations  of  oxygen  vacancy 
concentrations  on  the  Ca2'- 1-fold  coordinated  and 


T  °C 

U00‘C  1200T  1000‘C  800  X  600‘C 


Fif.  Sl  lof  fC  n.  1/T  plot 


Ca2‘-2-fold  coordinated  anion  sublattices;  V02  (1)  and 
Vo2  (2).  with  total  oxygen  vacancy  concentra¬ 
tion  [Vo2]  total  =  x  at  600°,  1000°,  and  1400°C 
calculated  according  to  Eq.  [6]  and  [8],  Also 
included  in  Fig.  7  are  the  site  fractions  of  1-  and 
2-state;  Nx  and  N2.  It  is  interesting  to  note  that  the 
variations  of  [Vo2  11  and  [Vo2  2]  with  x  is  steep¬ 
est  at  i  =  0.0625  and  become  more  gradual  with  in¬ 
creasing  x.  The  equiconcentration  point  [V02-  1]  = 
[Vo2  2]  shifts  to  the  higher  x  value  position  as  tem¬ 
perature  increases,  but  remains  in  the  low  x  region  in 
these  temperature  ranges. 

Experiment  shows  that  the  activation  energy  for 
ionic  conduction  increases  with  increasing  x  and  de¬ 
creases  with  increasing  temperature  at  constant  x.  Both 
features  are  accounted  for  by  the  present  model.  These 
follow  because  with  increasing  x,  the  relative  contribu¬ 
tions  from  a  and  a2.2  which  have  activation  en- 
laaS 

ergies  higher  than  a m  increase,  and,  as  shown  in  Fig. 
7  at  constant  x,  the  higher  the  temperature,  the  higher 
the  relative  contribution  from  [Vo2  1]  and  hence  the 
contribution  from  <ri-i,  which  has  the  lower  activation 
energy,  increases. 

Oxygen  diffusion  measurements  on  this  system  re¬ 
ported  so  far  are  those  by  Kingery  et  at.  (19),  Simp¬ 
son  and  Carter  (20),  and  Hagel  (21).  References  (19) 
and  (20)  report  agreement  of  the  measured  oxygen 
diffusion  coefficient  D0ia°  with  that  calculated  from 
ionic  conductivity  assuming  an  oxygen  transport  num¬ 
ber  of  1  and  the  correlation  factor  of  0.65  for  the  fluo¬ 
rite  structure.  In  contrast,  Hagel  (21)  obtained  an 
oxygen  diffusion  coefficient  D0is°  appreciably  lower 
than  that  calculated  from  his  own  conductivity  data. 
If  the  ionic  conduction  and  therefore  the  oxygen 
diffusion  in  this  system  are  multimode  processes  as 
described  by  the  present  model,  the  penetration  pro¬ 
file  of  O1**  from  the  gas  phase  would  not  yield  a  curve 
characterized  by  one  error  function  complement.  Such 
experiments  are  quite  difficult  to  carry  out  and.  de¬ 
pending  on  the  relative  contributions  of  each  mode,  the 
result  may  appear  to  fit  a  single  error  function  comple¬ 
ment  as  was  shown  for  the  diffusion  of  sulfur  in  NiO 
by  Howng  and  Wagner  (22). 
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Fig.  7.  Oxygen  vacancy  distribution  between  Ca2'-)  -fold  and 
Co^'-Z-fold  coordinated  anion  sublattica  as  a  function  of  total 
oxygen  vacancy  concentration  [Vo2]  total  =  x.  T  =  600°, 
1000°,  1400° C  [Vo2  1]:  oxygen  vacancy  concentration  on  Ca2'-1- 
fofd  coordinated  aniea  sablattice.  [Vo-  2]:  oxygen  vacancy  con¬ 
centration  ot  Ca2'-2-fold  coordinated  anion  sublattice.  Nt:  site 
fraction  ot  Co2'- 1 -fold  coordinated  anion  sublattice.  Nj:  site  frac- 
tien  of  Ca4'-2-fold  coordinated  anion  sublattice. 

Internal  friction  measurements  on  this  system 
conducted  by  Wachtman  and  Corwin  (.23)  seem  to 
support  the  above  speculation  on  the  oxygen  diffusion 
coefficient  and  partially  substantiate  our  proposed  theo¬ 
retical  model.  Wachtman  and  Corwin  observed  the 
symmetrical  peaks  associated  with  Vo2'  motion  in  the 
cubic  fluorite  phase  region,  but  the  peaks  were  too 
broad  to  result  from  a  relaxation  process  with  a  single 
relaxation  time.  They  thought  that  the  concentrations 
of  V02  and  Ca2-  are  so  high  that  many  complex  group¬ 
ings  and/or  various  arrangements  of  Ca2'  and  Vo2' 
occur  which  cause  the  various  types  of  local  V02 
jumps  and  therefore  a  distribution  of  relaxation  times. 
Their  interpretation  is  similar  to  that  of  Kroger  (.10) 
mentioned  previously. 

As  was  shown  in  the  foregoing,  the  proposed  theo¬ 
retical  model  seems  to  succeed  to  some  extent  in  giving 
a  quantitative  character  to  these  qualitative  interpre¬ 
tations.  The  basic  feature  of  the  theoretical  model  is 
that  it  extends  the  conventional  concept  of  dopant 
(impurity) -defect  complex  or  associate  for  highly  de¬ 
fective  systems  and  shows  in  a  straightforward  way  the 
occurrence  of  a  maximum  of  the  ionic  conductivity 
and  a  multimode  transport  mechanism  for  these  sys¬ 
tems.  However,  in  formulating  the  theoretical  model, 
crude  assumption  for  the  effective  anion  site  coordi¬ 
nation  number  of  Ca2'  to  the  second  nearest  neighbor 
<  Ca2'- V0-  )  configuration  was  made.  So  far,  our  knowl¬ 
edge  about  defect  interaction  at  the  second  (and  the 
higher  order)  nearest  neighbor  configuration  is  quite 
lacking.  Moreover,  the  statistical  fluctuations  of  the 
system  were  also  neglected.  (See  Eq.  [2],  Even  in  the 
composition  range  0.0625  ^  x  ^  0.125.  there  exist  also 
small  numbers  of  Ca2'-0-fold,  3-fold.  4-fold  anion 
sites,  which  were  all  neglected.  Also,  the  existence  of 


nonintegrally  Ca2'  coordinated  sites  such  as  Ca2'-l/3- 
fold,  Ca2'-5/3-fold  anion  sites  were  neglected.)  These 
fluctuations  would  be  dependent  on  the  Ca2'  distribu¬ 
tion  and  should  be  important  in  dealing  with  the  order- 
disorder  transition  in  this  system. 

Summary 

A  quantitative  theoretical  model  to  explain  the  vari¬ 
ations  of  ionic  conductivity  of  calcia-stabilized  zirconia 
with  oxygen  vacancy  concentration  and  temperature 
has  been  presented.  The  theory  extends  the  conven¬ 
tional  concept  of  the  impurity  (dopant)  -defect  complex 
or  associate  for  highly  defective  systems  and  shows 
that  ionic  conduction  in  this  system  proceeds  by  a 
multimode  mechanism:  Eq.  [13],  The  theory  predicts 
the  occurrence  of  a  rather  sharp  maximum  of  ionic 
conductivity  at  oxygen  vacancy  concentration  x  — 
0.0625  (12.5  m/o  CaO),  at  which  the  entire  anion  sub¬ 
lattice  is  just  covered  by  onefold  associated  (Ca2’-Vo2  ) 
complexes,  in  accordance  with  experimental  observa¬ 
tions. 

Values  of  characteristic  parameters  such  as  the  dif¬ 
fusion  coefficients  and  mobilities  of  oxygen  vacancies 
V02  through  and  between  1-fold  and  2-fold  Ca2'  co¬ 
ordinated  anion  sublattice;  Dm-d,  m h>,  and  the  equi¬ 
librium  constant  K  tor  the  distribution  ot  oxide  ions 
and  oxygen  vacancies  between  1-fold  and  2-fold  Ca2’ 
coordinated  anion  sublattice  were  obtained  by  a  curve- 
fitting  procedure  of  the  theoretical  expression  for 
ionic  conductivity  to  experimental  »(x,  T)  data  re¬ 
ported  in  the  literature.  Resultant  theoretical  curves  of 
isothermal  log  <j  vs.  oxygen  vacancy  concentration  x, 
inside  the  stability  region  of  the  cubic  fluorite  phase 
0.0625  x  =£  0.10  and  600°C  T°C  =£  1400°C  were 
shown  to  reproduce  satisfactorily  the  experimental 
curves  reported  by  various  investigators. 
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FIGURE  2 
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Parabolic  Plot  of  Oxidation  of  Nickel 
in  0.2  atm  O2  at  800°  C 

(TOP)  Schematic  diagram  of  location  of  P+  marker  after 
oxidation  in  0.2  atm  O2  at  800°C.  (BOTTOM)  The  value 

of  k  is  consistent  with  literature  data. 
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CORROSION  MECHANISM  IN  PURE  SO,  (2) 


•ni"1- *10 

NiSy(l) 

S02(9) 


-5 


REACTION  AX  SULFIDE/ OXIDE  Interface; 

S02  -  02  +  1/2S2 

p'  a:  3.0x10  ^  (atm)  For  Sulfide  Formation 

2  Nwo  V-l 

7Ni  +  S02(g)  -  4NI0  +  Ni-jSjCs  ,1) 

PgQ  x  0.1  (atm)  When  a^-0.12 

PgQ  a  5.3x10  ^  (atm)  When  a^-1,0 
2 

FIGURE  5.  Possible  mechanism  for  corrosion  of  nickel  in  S02_Ar  gas 

mixtures  at  800°C.  The  rate  determining  step  for  the  observed 
linear  kinetics  is  the  dissociation  of  S02(ads)  on  the  sulfide 
inner  layer.  L 
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FIGURE  6. 


Time,  (min) 

WEIGHT  GAIN  VS  TIME  FOR  NICKEL  AT  800°C 

Kinetics  for  207,  02"variable  S02"balance  argon  at  800°C.  Note 

the  two  regimes,  one  s  207o  in  which  the  corrosion  is  protective 
and  one  a  307,  in  which  the  kinetics  are  extremely  rapid. 
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WEIGHT  GAIN  VS  TIME  FOR  NICKEL  AT  800°C 


FIGURE  8.  Kinetics  in  307,  S02~variabLe  02_balance  Ar  at  800°C.  Note  the 
initial  trend  for  2,  5  and  107,  0„  and  the  behavior  for  s  207,  0 
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2h  Pre-ox 


Pre-oxidized  in  1  atm  0£  at  800  C 


•  Exposed  to  30%  SO2  -  70%  O2 
o  Exposed  to  30%  SO2  —  70%  Ar 


6h  Pre-ox 


6h  Pre-ox 


50 '' 200  250  300  350  400 

Time,  (min) 


Weight  gain  after  exposure  to  the  Cas  Mixtures  at  800  C 

FIGURE  9.  Effect  of  Pre-oxidation  in  pure  O2  on  corrosion  of  nickel  in 
S02*"02  ®aS  mixCures*  Samples  exposed  to  1  atm  O2  at  800  C. 
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FIGURE  10. 

Comparison  of  theory  (solid  lines)  with  experimental 
data  for  isothermal  ionic  conductivity  in  yttria- 
stabilized  zircon ia  as  a  function  of  vacancy  concen¬ 
tration. 


Experimental  Data: 

V  Dixon,  et  al.  (1963)  ref.  (6) 

0  Strickler  and  Carlson  (1964)  ref.  (7) 
X  Strickler  and  Carlson  (1965)  ref.  (8) 
□  Bauer le  and  Hrizo  (1969)  ref.  (9) 

A  Cassleton  (1970)  ref.  (10) 

•  Ioffe  et  al  (1978)  ref.  (11) 
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FIGURE  11. 

Typical  penetration  profile  for  sulfur  -  35  and 
tellurium  -  127m  in  calcia  stabilized  zirconia 
(15  m/o  CaO). 
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